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ABSTRACT 

^ A model  is  developed  for  use  in  sonar  performance  pre- 
diction for  systems  using  bistatic  echo  ranging.  A mathematical 
model  of  the  ocean,  suited  to  performance  calculations,  is  pre- 
sented. The  theoretical  pattern  for  a rectangular  array  is  de- 
veloped from  first  principles,  and  the  echo  level  and  background 
level  computations  are  outlined.  This  report  is  intended  as  a 
guide  for  future  performance  prediction  work  at  TRG.  t 
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SECTION  1 
INTRODUCTION 

This  report  covers  the  modeling  phase  of  work  done  at 
TRG  on  the  problem  of  performance  estimation  for  a sonar  system 
consisting  of  a physically  separate  source  and  receiver.  This 
work  is  in  the  nature  of  an  extension  of  the  previous  work  in 
performance  prediction  for  a coincident  source  and  receiver. 

The  goal  of  this  study  is  to  obtain  an  estimate  of  the  potential 
sonar  performance  of  systems  with  remote  receivers  for  applic- 
ation to  sonar  systems  design  and  application.  A situation  in 
which  a sonar  signal  is  received  by  a remote  receiver  can  be 
useful  enough  to  warrant  consideration. 

The  concern  of  this  study  of  sonar  performance  pre- 
diction in  the  bistatic  case  is  the  computation  of  the  echo-to- 
background  ratio  (E/B)  developed  by  a separated  source  and  re- 
ceiver system  against  a submarine  target.  Comparisons  of  the 
E/B  obtained  by  various  systems  provide  guidance  in  system 
design  and  application. 

After  a discussion  of  the  motivations  for  a bistatic 
eclio  ranging  study,  this  report  will  proceed  by  outlining  the 
ocean  model  used  in  the  systems  calevilations . Not,  the  echo- 
level  calculations  for  bistatic  situations  will  1 discussed  as 
an  extension  of  the  monostatic  echo- level  calculation  (as  al- 
ready Implemented  at  TRG) . Finally,  the  estimation  of  boundary 
reverberation  will  be  discussed  in  terms  of  a continuation  of 
the  techniques  now  u.sed  in  monostatic  reverberation.  This  re- 
port shows  the  direction  in  which  sonar  performance  prediction 
for  separated  source  and  receiver  at  TRG  is  proceeding.  (The 
monostatic  calculations  are  set  forth  in  023-TM-66-16,  a com- 
panion volume  (Reference  24  of  this  report)). 

MOTIVATION  FOR  THE  BISTATIC-ECHO  RANGING  STUDY 


The  potential  advantage  of  a sonar  system  with  a re- 
mote receiver  over  a system  with  coincident  source  and  receiver 


win  1-2-4 100 ‘J 


is  due  to  its  added  freedom  in  the  placement  of  the  receiver.  This 
advantage  is  due  not  only  to  the  potentially  higher  echo-to-back- 
ground  ratios  obtainable,  but  also  to  the  geometric  situation 
which  may  allow  a more  precise  determination  of  the  target  lo- 
cation. A sonar  system  restricted  to  monostatic  echo  ranging 
must  be  content  with  whatever  propagation  conditions  exist  be- 
tween it  and  the  target  for  both  the  transmission  and  reception 
paths.  On  the  other  hand,  a bistatic  system  may  take  advantage 
of  more  favorable  propagation  conditions  for  the  target- to-re- 
ceiver  path.  When  the  receiver  is  positioned  so  that  the  re- 
ception path  is  direct,  the  target  localization  is  essentially 
independent  of  the  bottom  rou;_,hness,  which  presents  a problem 
for  bottom  bounce  sonars.  (In  addition,  multiple  receivers  may 
be  used  to  pin-point  the  target.) 

Furthermore,  a system  with  a remote  receiver  may  lead 
to  more  advantageous  target  localization  than  monostatic  systems. 
With  the  position  of  the  source  and  receiver  known,  the  echo 
travel  time  determines  a contour  on  which  the  target  must  lie. 

(The  specification  of  this  contour  is  discussed  in  Section  VI.) 

To  fix  the  target  location  uniquely,  further  information  is 
needed.  A bearing  from  the  source  (which  also  receives)  could 
locate  the  target . 

Alternatively,  the  receiver  covild  also  contain  a 
transponder  to  take  advantage  of  the  Loreli  target  localization. 
(See  Reference  30.)  A Loreli-type  system,  with  a remote  re- 
ceiver could  localize  the  target  as  the  intersection  of  the 
Loreli  bearing  with  the  contour  determined  by  the  echo  travel 
time.  (The  Loreli  bearing  is  found  by  subtracting  the  sonar 
vector  from  the  source  to  the  receiver/transponder  from  the 
sonar  vector  to  the  target,  which  is  also  detected  monostatically .) 
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This  subtraction  serves  to  eliminate  systematic  errors  (e.g. 
velocity  gradients  in  the  horizontal  direction,  or  sonar  align- 
ment errors)  from  the  target  localization.  As  with  Lorell 
systems,  a sonar  system  with  a remote  receiver  in  the  vicinity 
of  the  target  converts  the  long-range  fire  control  problem,  re- 
quiring high  absolute  accuracy,  to  a short-range  fire  control 
problem,  requiring  only  relative  accuracy.  The  geometry  is 
illustrated  below  in  Figure  1-1. 

Several  studies  have  indicated  the  superiority  of 
multistatic  echo-ranging  systems.  (See,  for  example.  Reference 
31) . Work  in  underwater  target  classification  has  shown  that 
bistatic  techniques  produce  more  infoi*mation  than  monostatic 
systems  (because  the  echo  structure  depends  on  the  bistatic 
angle  (defined  in  Section  V)) . Current  work  in  target  class- 
ification is  based  on  the  use  of  multistatic  systems  which  can 
provide  the  necessary  information.  Theoretical  studies  were 
verified  by  at-sea  experiments  in  locating  and  classifying  sub- 
marine targets,  lndica*’tng  the  f•'?»^^billty  of  multistatic  sys- 
tems in  target  classi f ication . Furthermore,  the  tactics  re- 
quired for  a multistatic  classification  position  the  elements 
of  the  ASW  system  so  that  they  can  maintain  contact  with  a fast, 
maneuverable  target. 

In  addition,  studies  and  experience  show  that  a multi- 
static system,  because  of  its  higher  information  rate,  can 
classify  a target  in  a shorter  time  than  a monostatic  system, 
although  the  multistatic  geometry  takes  longer  to  set  up  . The 
multistatic  system  has  further  advantages.  First,  it  receives 
more  information  on  the  target's  speed  and  course  due  to  the 
additional  .sensors  in  the  geometry.  Second,  the  tactical  dis- 
bursement of  the  ASW  system  elements  restricts  the  target  es- 
cape capabilities. 

In  addition  to  studies  based  on  multistatic  echo 
ranging,  blstatic  situations  are  commonly  encountered  in  cur- 
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Figure 


TARGET 


1-1.  BISTATIC  TARGET  LOCATION  GEOMETRY 
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rent  tectlcel  situations.  Convoying  operations  offer  an  illus- 
tration of  the  effectiveness  of  a remote  receiver  system. 

Here,  a number  of  escort  vessels  operate  in  the  active  mode. 
Monostatic  operation  restricts  each  ship  to  listen  only  to  Its 
own  echo.  An  understanding  of  the  blstatic  problem  would  allow 
escort  vessels  to  take  advantage  of  the  other's  transmissions. 
The  use  of  sonobuoys  to  receive  active  transmissions  (from  some 
other  source)  also  Involves  a blstatic  situation. 

The  above  examples  Indicate  the  possible  advantages 
of  blstatic  sonar  systems  Implying  a need  for  a performance  pre- 
diction capability  for  blstatic  sonars.  Such  a capability 
would  provide  guidance  In  system  design  and  application.  Fam- 
iliarity with  the  blstatic  problem  would  not  only  allow  more 
tactical  freedom  In  the  use  of  sonars,  but  it  wovild  also  shed 
light  on  the  operational  problem  of  mutual  sonar  interference. 
Optimal  tactical  positioning  of  sonar  ships  requires  an  inves- 
tigation of  the  Interference  problem,  Involving  the  same  cal- 
culations as  the  blstatic  echo-ranging  problem. 

For  system  performance  purposes,  the  general  problem 
concerns  the  echo-to-background  ratio  Cor  arbitrary  placement 
of  the  source  and  receiver.  The  oxtenslvm  of  present  prediction 
capabilities  to  the  blstatic  case  will  allow  an  Investigation  of 
the  effects  of  ship  motion,  incUuling  translation,  roll,  pitch, 
and  yaw,  on  sonar  performance,  as  a special  case.  (Different 
transmit  and  receive  systems  on  the  same  ship  could  also  be 
evaluated^  Furthermore,  the  extension  of  the  blstatic  solution 
to  the  multlstatlc  problem  for  sonobuoy  arrays  and  multiple  es- 
cort vessels  configurations  Is  straightforward.  For  the  reasons 
outlined  above,  an  investigation  of  blstatic  echo  ranging, 
from  the  point  of  view  of  performance  prediction,  is  needed  for 
a more  complete  understanding  of  the  operations  of  complex, 
multi-element  sonar  system. 

An  understanding  of  the  various  bLstatic  sonar  ap- 
plications described  above  should  lead  to  Improved  utilization 
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of  current  and  planned  sonar  systems.  In  addition,  a remote 
passive  receiver  will  not  alert  a target  to  the  fact  that  it 
has  been  detected;  an  unalerted  target,  coupled  vlth  improved 
localization  should  result  in  greatly  improved  weapon  effective- 
ness. Finally,  as  a harassing  tactic,  the  possible  use  of 
remote  receivers  would  be  an  added  threat  for  the  submarine 
commander  by  restricting  his  use  of  good  propagation  paths, 
degrading  his  infornution  and^  hence,  his  performance. 
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SECTION  II 
GENERAL  CONCEPTS 


The  fundamental  concepts  used  In  systems  analysis  are 
discussed  In  Section  I of  Reference  I.  Some  of  the  topics  cov- 
ered In  this  section  are  summarized  below; 

Directivity  Index.  PI; 

Att  X Effective  Area 

2 

(Wavelength) 


1 yd)  =s  101.6  + 10  logj^Q[ Power 

Out(kw) ]+DI. 

Noise  Level, 

If  the  noise  is  truly  isotropic,  the  eauivalent  plane 
wavo  noise  level,  given  by; 

Lepy(db  re  lubar  ) = Isotropic  spectrum  level  (db  re 

lubar^/cps)  - DI  f 10  log^O 
[ bandwidth (cps) ] 

★ 

If  the  noise  is  uncorrelated  between  elements,  e.g.,  flow  in 
duced  noise,  then  L^p^  is  given  by: 

2 

^epw  “ Single  element  spectrum  level  (db  re  li^ibar  /cps) 
-10  logj^Q(t\o.  of  elements)  + 10  log j^q (bandwidth) 

Background  Level,  L^ 

The  background  level  is  assumed  to  be  cocposed  of 
reverberation  and  noise.  Since  these  background  signals  are  un- 
correlated, the  total  power  is  the  sum  of  the  powers  of  the 
signals . 


Dl(db)  = 10  logj^o 

Source  Level,  L^ 

L (db  re  1 ubar  at 


Over  the  range  of  time  delays  used  for  beamforming 
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Echo  Ltfvoi.  Lg 

^ocho  Uva.  «re  computea  fron>  the  etsndora  echo 
ranging  equation, 

L^(ab  ce  laber^)  - ' (trenecleelon  t 

reception  deviation  Iobb) . 


where;  N„  - target  eerength  (In  db) . 

g - two-way  propagation  lose  (In  db) . 
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SECTION  III 
OCEAN  MODEL 


The  basic  ocean  model  assumed  for  the  echo  and  rever- 
beration level  calculations  Is  described  In  the  following  sec- 
tions. 

A . PROPAGATION 

1.  Propagation  Loss 

The  propagation  loss  is  given  by  the  sum  of  the  In- 
dividual ray  path  losses;  namely,  propagation  loss  = spreading 
loss  + reflection  losses  + absorption  loss. 

The  spreading  loss  In  an  Isoveloclty  medium  given 
by  20  (slant  range) . The  calculation  of  spreading  loss 

In  a refractive  medium  Is  described  below.  The  reflection 
processes  are  discussed  below. 

2.  Absorption  Loss 

The  absorption  loss  coefficient  (db/kyd)  Is  repre- 
sented by  a general  power  law  equation;  that  Is,  the  coeffi- 
cient Is  equal  to  a constant  times  the  frequency  raised  to  some 
power.  The  absorption  loss  in  db  is  given  by  the  product  of  the 
coefficient  and  the  slant  range  to  the  target  (which  is  approxi- 
mated by  the  mean  sound  velocity  x travel  time  along  the  path). 

3 . Ray  Tracing 

The  velocity  profile  is  assumed  to  be  a function  of 
depth  only,  and  Is  approximated  by  a series  of  layers,  each 
with  a constant  velocity  gradient.  For  a ray  traversing  a 
single  layer,  the  basic  formulas  used  for  the  ray  tracing  calcu- 
lations are  summarized  below; 


\R  = (sin  0.  - sin  0J 
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dAR 

39;' 


g 


I 


I 

sin  6^ 


(2) 


AT 


I (Usineg)  (l-sinO^) 

.■g  ^°8e  TTriTnS^TTT+sTnB^ 


(3) 


dAT  I d\R 
O X o 


(4) 


where : 


AR  Is  the  incremental  horizontal  range 

0 Is  the  Inclination  angle  of  the  ray  at  the 
° source  depth, 

0 is  the  inclination  angle  of  the  ray  at  its  entry 
® into  the  layer, 

0^  is  the  Inclination  angle  of  the  ray  when  it  leaves 


the  layer, 


at  is  the  incremental  travel  time  (in  seconds)  through 
the  layer, 

is  the  Snell's  law  constant  and  is  given  by 


o 

cosS, 


where : 


Vq  is  the  sound  velocity  at  the  source, 
g is  the  velocity  gradient  in  the  layer  and  is  given  by: 


g = ■ , where: 

V is  the  velocity  of  sound  at  the  depth  where  the  ray 
leaves  the  layer,  and 

Ad  is  the  thickness  of  the  layer. 


Some  of  these  quantities  are  illustrated  in  Figure  3-1. 

A suOTnatlon  over  the  layers  involved  produces  the  hori- 
zontal range  and  travel  time  for  a ray  path.  In  a non-refractive 
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meJiuin,  the  Intensity  at  a point  far  from  the  source  obeys  the 
Inverse  square  law:  I = analogue  of  the  Inverse 

square  law  for  the  refractive  medium  is  the  spreading  loss 
factor.  For  a ray  bundle  travelling  through  a series  of  re- 
fractive layers,  it  is  given  by: 


cos9^ 

— rrm 

3B^  “1"  Sp 


where ; 


k is  defined  as  the  ratio  of  the  sound  intensity  at 
a point  p to  the  source  intensity  referenced  to 
one  yard, 

R is  the  total  horizontal  range  (r  >'\Ri), 

1 ^ 

0Q  is  the  depression  angle  of  the  ray  at  the  source,  and 
9p  Is  the  depression  angle  of  the  ray  at  point  p. 

The  reader  Is  directed  to  references  2 through  A for  the  deriva- 
tions and  explanations  of  the  above  relationships. 

4.  Bottom  Loss  (Specular  Reflections) 


There  are  currently  two  seml-empir leal  equations  used 
for  the  calculation  of  the  specular  reflection  coefficient.  The 
first  one  Is  based  on  a theoretical  two-layer  model  for  an  absorb- 
ing bottom  as  reported  by  Mackenzie  in  Reference  3.  The  formula 
has  been  modified  by  the  addition  of  some  curve  fitting  parameters, 
enabling  a satisfactory  fit  to  the  modified  AMOS  bottom  loss  data 
(Reference  6) . For  oceanographic  areas  which  are  adequately  des- 
cribed by  the  two-layer  model,  the  free  parameters  may  be  suppressed, 
leaving  the  theoretical  expression. 

The  bottom  loss  equation,  for  the  specular  reflection 
coefficient,  r,  is 

V f.  1 CjU-CjO)  (») 

(T+PQ  Sin  «)^+  S-i  J 
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where: 


! 

, s 


I i 


- 

S2- 

A - 

B 
b 


(A-B)/2, 

(A+B)/2, 

(B2+4b2)^/2^ 

p2  cos^g-l+b^,  and 


Here, 

P le  the  longitudinal  velocity  ratio  (botton/water) , 

Q le  the  density  ratio  (botton/water), 

b is  related  to  the  absorption  coefficient  of  the 
bottora, 

f Is  the  frequency  In  cps,  and 

^2*  ^3*  ® curve  fitting  constants. 

For  example,  the  constants  used  to  fit  the  AMOS  data  are: 


P - 1.025 
Q - 1.2634 
b^  - 0.0381 
e - 0.0 


B 3500  cps 
C2  = 1.315 
- -5.04 


An  alternative  bottom  loss  expression  is  due  to  Dr. 
V.  Watson  of  NEL: 


Ng  -[6  + 22(p-0.27)  + 10  log 


tanh 


N 


373* 


\ ( 0 , 

"irrr 


» ^kHz  - ^ 


(7) 


where: 


Ng  is  expressed  In  db. 

p is  the  porosity, 

^kHz  frequency  In  kHz  and 

0 Is  the  grazing  angle  of  Incidence  In  degrees. 
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5 . Surface  Loss  (Specular  Reflection) 

There  is  insufficient  information  in  the  current  lit- 
erature to  allow  a reliable  estimate  of  the  surface  reflection 
coefficient  at  the  frequencies  used  in  bottom-bounce  sonars. 
References  26,  7,  R and  9 indicate  that,  because  of  the  large 
impedance  mismatch  between  water  and  air,  the  surface  reflection 
coefficient  is  close  to  unity.  Although  the  contributions  of 
these  effects  are  not  well-known,  the  surface  can  affect  the 
signal  through  the  following  mechanisms: 

1.  It  can  scatter  energy  out  of  the  ray  bundle.  (Reference 
26  indicates  that  energy  may  leave  the  surface  channel  by  re- 
flection from  wave  facets)  . 

2.  It  can  absorb  energy  irom  the  signal  (through  the  action 
of  near-surface  trapped  air  bubbles) . 

3.  It  can  degrade  the  phase  coherence  along  the  wavefront, 
since  the  signal  is  reflected  from  randomly  oriented  facets  of 
the  surface.  For  a highly  directional  sonar,  the  full  gain  of 
the  sonar  could  be  lost  for  a signal  propagating  over  a surface 
reflected  path  (as  an  upper  limit)  . This  cohert'nce  loss  applies 
only  to  the  echo;  reverberation  Is  generally  assumed  to  be  an 
incoherent  process,  so  the  reverberation  intensity  is  not  changed 
by  a loss  in  signal  coherence. 

When  estimating  sonar  performance  in  the  bottom  bounce 
mode,  only  a single  path  to  the  target  is  cvirrently  considered. 
This  should  provide  a conservative  estimate  of  the  sonar's  poten- 
tial performance  (since  it  implies  a surface  reflection  coeffic- 
ient equal  to  zero. 

B.  SCATTERING 

1.  Bottom  Scattering  (Non-specular  reflections^ 

The  bottom  scattering  coefficient  is  assumed  to 

obey  a generalized  Lambert's  Law  expression: 

m n 

^b  “ ‘■^bo  ®i  '• 

*See  also  References  34  - 37. 
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where 

is  a constant  of  the  diffusing  surface, 

cos  is  the  incident  direction  cosine  with  the  normal 
to  the  aurfacn, 

cos  0s  is  the  scattering  direction  cosine  with  the  sur- 
face  normal, 

> is  the  azimuthal  angle  between  the  incident  and 
scattering  directions,  and  m & n are  curve-fitting 
constants,  chosen  so  that  the  expression  agrees 
with  experimental  results.  (Isotropic  scattering 
corresponds  to  ro=n--0) 

2.  Sea  Surface  Scattering 

Very  few  experiments  have  been  reported  dealing  with 
the  spatial  dependence  of  non-specular  reflection  from  the  sea 
surface*.  At  present,  a cos  > azimuthal  dependence  will  be 
assumed,  as  in  bottom  scattering,  to  represent  experimental 
data  when  it  becomes  available.  An  extension  of  the  results  for 
reciprocal  ray  paths  is  noted  below.  Modifications  will  be  made 
as  further  experimental  evidence  becomes  available.  The  equation 
for  the  sea  surface  scattering  coefficient  ii  are  based  on  the 
experimental  work  reported  by  Chapman,  et  al,  in  Rr  _rences  7 & 8 
for  the  intermediate  and  large  grazing  angles  over  reciprocal  ray 
paths.  The  equations  are  given  below. 


(8) 


where  is  based  on  Eckart's  work  and  is  given  by 


se 


(9) 


where 


a*"  - 0.003  0.00263W 


and  u is  based  on  an  "eyeball"  fit  to  the  corrected*  Chapman 
and  Harris  data  and  it  is  given  by 

See  Reference  37. 
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u 


sc 


-(2.7  log^Qp+1.16)  g 0.33p 

10 


(10) 


where 


p - 158 


In  the  above  equations, 

W is  the  wind  speed  In  knots, 

9 Is  some  combination  of  the  incident  and  scattered 
grazing  angles 

(6^  and  9^,  respectively)  chosen  to  represent  experi- 
mental data  |e.g.  9 = ^ g — - ; or  9 = J,  and 

£ Is  the  frequency  in  Hertz. 

3 . Deep  Scattering  Laver 

Volume  reverberation  at  the  frequencies  used  In  bottom 
bounce  sonar  is  due  primarily  to  fish  with  resonant  or  near-reson- 
ant air  bladders.  For  any  given  location  and  time,  these  scat- 
terers  are  generally  confined  to  a relatively  narrow  depth  In- 
terval. They  may  be  treated  as  an  equivalent  scattering  layer, 
commonly  referred  to  as  the  deep  scattering  layer  (DSL) . The 
statistics  of  the  low  frequency  DSL  are  currently  under  Inves- 
tigation by  various  Navy  laboratories,  but,  at  present,  are  not 
well  known.  Qualitative  discussions  of  the  DSL  appear  in  Ref- 
erences 11  to  16;  most  of  the  quantitative  data  for  the  bottom 
bounce  frequency  range  is  contained  In  References  7,  17,  18,  and  38. 
Further,  reverberation  from  the  deep  scattering  layer  Is  due 
to  scattering  from  biological  material  with  effective  dimensions 
less  than  a wavelength  and  Is  generally  considered  to  be  Iso- 
tropic . 


B WP11-2-A1005 


SECTION  IV 

CONFORMAL/ PLANAR  ARRAY  CONSIDERATIONS 


A.  THE  PATTERN  FUNCTION 

The  exact  expressions  for  the  pressure  pattern  for  a 
rectangular  array  of  uniformly  driven  point  sources  is  derived 
from  basic  principles.  Then  adjustments  are  made  to  account  for 
the  effects  of  Individual  element  directivity  and  for  the  effect 
of  a finite  baffle, 


In  the  most  general  case,  the  pattern  function  (V)  is 

given  by; 


V-  f dv„ 

J Ro 

all  space 


where:  k •>  2rT/x  (1^  a unit  vector  in  the  observation  direction), 

ko"  2n/\  (1^^  a unit  vector  in  the  pointing  direction) 

Rq  is  the  vector  to  the  source  point,  and  p (R^)  is  the  source 
velocity  distribution.  Since  the  pattern  is  of  interest  only  in 
the  far  field,  the  vector  R,  is  written  R - R ♦ r.  In  the  far 
field,  R » r,  so  that  r need  be  retained  only  as  a phase  factor, 
allowing 


V 


'\J 


dV 


all  space 


i(k-kp)  • r 


When  all  the  transducers  are  driven  with  uniform  velocity, 
the  source  velocity  distribution  r(R„)  - L 8(Rq-Rj),  in  the 

•pproxloMtlon  where  the  transducers  may  be  considered  as  point 
sources.  The  voluow  Integration  over  the  Dirac  functionals  reduces 
the  pattern  function  to  a sum  of  exponentials. 


Starting  with  the  simple  case  of  a line  array  of  n 
uniformly  driven  point  sources,  spaced  a distance!  d, apart,  the 


HKy 


'Mr  w-af.V:-i».?ri5HEr-sewc:a4ta.-asa-rri!<teK.« 
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I 

I 

R 

r 


pattern  V is  given  by; 


V - V y c 
n»»l 


ini'. 


where  z •«  kd(co8^  - cos^/q),  and 

k 2ir  /^  is  the  wave  number. 


For  a line  array « the  array  axis  forms  the  symmetry 
axis  of  the  pattern,  ip  is  the  angle  formed  by  the  direction  of 
interest  with  the  line  of  the  array  (as  Illustrated  In  Figure  4-1) . 
Because  of  the  symmetry  of  the  pattern^  f may  be  specified  by  the 
half-angle  of  the  cone  generated  by  the  spatial  direction  and 
the  array  axis.  Alternatively » with  a view  toward  generalizing 
the  definition  of  lor  a rectangular  array > cos^'  Is  the  direc- 
tion cosine  of  the  spatial  direction  with  respect  to  the  line  of 
the  array. 


The  pattern  Is  a geometric  scries 


oe 


iz(l-e^!^) 


1-e 


iz 


''o*  ‘ " } \,izn:-i:zny 


which  Is  now  a familiar  form.  We  have,  finally, 

V » ^""^^^sln  nz/2 

sin  z/2 

which  Is  the  unnormallzed  pattern  for  a line  of  point  sources. 


B.  SHADING  THE  ARRAY 

The  above  discussion  holds  only  if  the  elements  of  the 
array  are  driven  with  uniform  velocity.  The  shape  of  the  pattern, 
may  be  altered  by  driving  the  array  non-uniform  velocity  distrib- 
ution. A common  application  Involves  controlling  the  velocity 
distribution  ("shading  the  array")  in  a special  fashion  to  reduce 
the  side-lobe  level,  at  the  expense  of  broadening  the  main  beam. 

To  model  this  (or  other)  velocity  distribution  without  having  to 
evaluate  the  sunmiatlon  for  V explicitly,  the  array  is  considered 


i 
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to  be  co»Tiposite,  as  shown  In  Figure  4-2. 

In  Figure  4-2,  Che  bj's  measure  Che  displacement  of  the 
sub-arrays  from  the  sub-array  Chat  represents  the  uniform  velocity 
distribution  (subscript  "o").  The  amplitude  of  each  sub-array  is 
denoted  by  the  Vj . The  velocity  distribution  is  approximated  by 
superimposing  several  sub-arrays  so  that  the  velocities  of  the 
sub-arrays,  at  any  point  along  the  array,  add  to  the  desired  vel- 
ocity . 

It  follows  that  the  pattern  from  the  shaded  (composite) 
is  given  by  the  superposition  of  the  patterns  from  the  sub-arrays 

V - T V ”45''/2  ikb  z/d 

^ slh  z-h  ® ^ • 


Consequently,  Che  unnormuiized  pattern  function  for  a shaded  line 
array,  composed  of  J sub-arraya,  is  given  by 


lz/2 

e 

8inz/2 


J i2(n  /2  + K, ) 
y Vjg  J ^ 8ln(njZ/2) 


where  the  J sub-array  starts  with  the  element  Kj  (numbered  from 
the  left  end) , 


Normalization  is  afforded  by: 


Vjj  - V (z-0) 


J5o 


C.  PATTERN  FOR  A RECTANGULAR  ARRvW 

The  above  discussion  has  been  for  a line  array.  The  ex- 
tension of  these  results  to  a rectangular  array  makes  use  of  the 
well-known  Product  Theorem,  The  pattern  of  a rectangular  array 
of  uniformly  driven,  equi-spaced  elements  may  be  obtained  by  multi 
plying  the  pattern  of  the  line  array  composed  of  one  row  of  the 
elements  of  the  rectangular  array  by  the  pattern  of  the  line  array 
composed  of  one  column  of  the  elements. 
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OHt«IN  rOR  CLCMCNT  ENUMIRATION  DISTANCE  ALONG  ARRAY 


FIGURE  ^♦-2.  APPROXIMATE  VELOCITY  DISTRIBUTION  FOR  A SHADED 
LINE  ARRAY 
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The  Fioduct  Theorem  is  based  on 


Ikr 


The  height  and  length  patterns  (of  the  line  arrays  re- 
ferred to  above)  are  formed  separately.  Ignoring  the  phase 
angle  which  lacks  physical  significance,  the  pattern  of  a uni- 
formly driven,  uniformly-spaced  rectangular  array  (the  length 
I and  height  spacing  may  be  different)  is  given  by  V = Vj^.  Vjj, 

I 

; with  and  computed  as  described  above. 


D.  COMPUTING  THE  DIRECTION  COSINES 

The  preceding  discussion  used  direction  cosines  with 
respect  to  the  array  to  define  the  cosines  appearing  In  the 
expression  z kd(cos'^'  - cos^/^)  . The  arguments  for  the  length 
and  height  patterns  are  defined  In  terms  of  a coordinate  system 
relative  to  the  array.  In  practice,  one  is  concerned  with  sonar 
coordinates,  with  depression  angle  measured  positively  above  the 
horizontal,  and  azimuth  measured  from  the  ship's  heading.  This 
is  illustrated  in  Figure  4-3. 

The  direction  cosines  relative  to  the  array  are  ob- 
tained from  the  spatial  direction  cosines  (defined  relative  to 
sonar  coordinates)  by  a rotation  of  coordinates  using  the  Euler 
angles  defined  as  follows:  The  yaw  angle  Y is  defined  In  the  same 
manner  as  the  azimuth  angle.  The  roll  angle  R is  taken  as  de- 
noting rotation  about  the  ship's  longitudinal  axis,  and  the  pitch 
angle  P,  as  measuring  rotation  about  the  ship's  transverse  axis. 

The  rotation  of  coordinates  is  accomplished  by  a m^itrix 
multiplication;  z A'C_, 

where  c is  the  vector  of  spatial  direction  cosines. 


(cos  0 cos  0 
cos  0 sin  0 
sin  0 
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Combining  the  rotations. 


cosP  cosY-cosR  sinY  slnp 
-stnP  cosY-cosR  slnY  cosP 
i sinR  sinY 


cosP  sinY+cosR  sinY  cosP 
“SinP  sinY+cosR  cosY  cosP 
-sinR  cosY 


slnP  sinR 
cosP  sinR 
cosR 


where 

R«  roll  angle 
Y«  yaw  angle 
P=  pitch  angle 


As  an  illustration,  the  calculation  of  the  pressure 
pattern  for  an  array  with  no  pitch  and  yaw,  with  a roll  angle 
R a -0^  is  shown  below.  (This  is  the  pattern  function  presently 
used  at  TRG  for  the  C/P  array  calculations). 

Expanding  z - we  have: 


cosO  COS0 

cos0,j,cos0stn0-stn0^8in6 

sln0,^cos08in0+cos0,^8in6, 


1 

0 

0 


0 0 
cos0,j,-sin0^ 
sin0,^  cosQji 


COS0COS0 

cosOsin0 

sin9 


Taking  the  first  and 

as  and 


third  components  of 


the  c vector 


SinL  ?\  (cOSO  CO80-COs9n  COS0ri). 
ain[— — (cos0  cos0-cos9_  cos0  ) 

• A O O 


sin 

^nRirdn  j 

|co89.j.(sin9-sin9^)+sin9,j,(cos9sin0-cos9^sin0^) 

A ' 

sin 

1 

|(cos9.j,(sin9-sin0^)+sin0^(cos9sin0-cos0^sin0^) 

11 

1 1 

L A ' 

/J 

after  introducing  the  pointing  direction  (zero,  subscript) . This 
expression  is  consistent  with  th'i  expression  given  in  Reference  24. 


4-8. 


WPLl-2-41005 


E.  INDIVIDUAL  ELEMENT  DIRECTIVITY  FACTOR 

The  above  calculations  make  the  assumption  that  the 
C/P  array  is  composed  of  point  sources.  The  directivity  of  the 
pistons  of  finite  diameter  may  be  taken  into  account  by  the 
factor  Q = (ka  sin  (w)  ) /ka  sin  (w) 

where  Jj^  denotes  the  first  order  cylindrical  Bessel 

function, 

k is  the  wave  number,  as  before, 
a is  the  piston  radius,  and 
w « cos  ^(22) 

(22  is  the  second  component  of  the  z vector  introduced  above,  the 
direction  cosine  with  respect  to  the  normal  of  the  array).. 

Noting  that  ■hdi-  U^(t),  the  zero  order  associ- 

ated Tchobycheff  polynomial,  the  individual  element  directivity 
factor  can  be  written  in  the  computationally  convenient  form:  Q =- 

4.vr?)/  ka  ■Clf’  (noting  Uy(t)  = 1) 

^For  square  elements,  a similar  expression  may  be  used), 

F.  FINITE  BAFFLE  EFFECT 

The  above  expressions  ai-o  Lor  on  array  In  a free  acoustic 
field.  As  mounted  on  a ship,  the  transducers  are  in  a rigid 
baffle  (after  resonant  stiffening)  which  has  the  effect  of  doub- 
ling the  pressure  pattern  (since  radiation  is  restricted  to  a 
half-space).  Diffraction  is  a further  effect  of  the  finite  baffle. 
The  form  of  the  expression  for  the  baffle  diffraction  has  been 
reported  in  several  IRG  reports  under  this  contract.  Reference 
28  developed  equations,  in  terms  of  the  Fresnel  integrals,  for 
the  finite  baffle  effect  for  baffles  large  with  respect  to  a wave- 
length. Reference  29  compared  the  patterns  of  an  array  on  a cylin- 
der with  the  pattern  of  an  array  on  a plane  and  found  that  the 
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effects  of  curvature  are  small.  A computationally  convenient 
form  of  the  finite  baffle  diffraction  expression  Is  a subject 
for  further  Investigation. 

G.  ARRAY  SELF-NOISE 

Derived  from  at-sea  measurements  (see  Reference  20), 
the  noise  level  used  In  systems  calculations  Includes  flow  noise, 
machinery  noise,  vibration,  and  ambient  sea  noise.  Adding  ele- 
ment outputs  In  an  array  Increases  noise  and  signal  at  different 
rates  due  to  an  averaging  of  the  random  noise  over  the  many 
elements  of  the  C/P  array  while  the  non-stochastlc  echo  adds 
coherently.  Accordingly,  the  array  gain  against  noise  Is  pro- 
portional to  the  number  of  elements  in  the  array. 

Array  Gain  ^ 10  (number  of  elements  in  array)  . 
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SECTION  V 

DETEEIMINATION  OF  THE  ECHO  LEVEL 

The  echo  level  is  the  acoustic  intensity  at  the  re- 
ceiver due  to  the  reflection  of  energy  from  the  target.  The 
basic  process  concerns  the  flow  of  energy  from  a point  "A"  in 
a refractive  medium  to  a point  "B".  The  bistatic  echo  is  de- 
veloped in  two  stages.  The  energy  produced  at  the  source  is 
transmitted  over  some  ray  path  to  the  target.  Then  a fraction 
of  this  energy,  given  by  the  target  strength,  is  directed  over 
another  ray  path  to  the  receiver.  A knowledge  of  these  ray 
paths  permits  the  calculation  of  the  propagation  loss,  the  de- 
viation losses  due  to  Che  array  patterns,  and  the  insonification, 
all  of  which  are  considered  below.  Finally,  these  results  are 
combined  with  the  target  strength  Co  yield  the  echo  level. 

In  general,  several  types  of  propagation  paths  may 
exist  between  the  source  and  the  target,  and  between  the  target 
and  receiver.  A given  transmission  and  reception  path  pair  will 
create  an  echo  at  the  time  t'.  Other  combinations  of  paths  will 
deliver  the  echo  at  some  different  time  t'',  and  may  be  consid- 
ered separately.  To  simplify  the  following  discussion,  results 
are  presented  for  any  particular  pair  of  ray  paths.  Since  the 
nature  of  the  echo  level  is  independent  of  the  details  of  Che 
ray  paths,  the  discussion  holds  for  any  combination  of  trans- 
mission and  reception  ray  paths.  If  multiple  paths  exist  such 
that  the  respective  echoes  overlap  in  time,  then  the  observed 
echo  is  a summation  of  each  of  the  contributing  echoes.  All 
processes  considered  below  have  this  implied  summation  over  all 
transmission- reception  path  pairs.  (This  results  in  a consid- 
erable simplification  of  the  equations;  supplying  the  summation 
over  all  ray  paths  is  a simple  operation.) 
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A.  RAY  TRACING  (FINDING  THE  RELEVANT  RAY  PATHS) 

This  section  is  concerned  with  the  determination  o£ 
the  source  depression  angles  for  the  ray  path  starting  at  the 
point  "A”  and  ending  at  the  point  "B**  in  a refractive  medium. 

The  transmission  and  reception  paths,  specified  by  their  de- 
pression angles,  are  found  as  described  below. 

Section  III  spells  out  the  ray-tracing  equation  for 
an  ocean  modelled  by  a refractive  medium  consisting  of  horizon- 
tal layers  with  a constant  velocity  gradient.  For  the  purposes 
of  systems  calculations,  the  sound  velocity  is  adequately  re- 
presented as  a function  only  of  depth.  While  the  details  of  the 
results  below  depend  on  this  representation  of  the  ocean,  the 
conclusions  and  techniques  are  valid  for  any  convenient  repres- 
entation of  the  velocity  profile. 

The  technique  of  ray-tracing  discussed  in  Section  III 
allowed  the  calculation  of  the  horizontal  range  traversed,  given 
the  source  depression  angle.  The  echo  level  computation  requires 
the  solution  of  the  inverse  problem:  Given  the  horizontal  range 
to  the  target,  find  the  source  angle  of  the  ray  that  reaches  the 
target.  Since  the  formula  for  the  horizontal  range  travelled  by 
a ray,  as  a function  of  its  source  angle,  does  not  allow  inver- 
sion to  yield  an  analytic  expression,  iterative  techniques  are 
used  to  find  the  source  depression  angle  for  the  ray  connecting 
points  "A”  to  "B". 

By  their  nature,  iterative  schemes  require  an  initial 
estimate.  For  bottom-bounce  ray  paths,  source  angle  estlmiites 
are  derived  from  Isovelocity  considerations.  Methods  for  de- 
veloping estimates  of  source  angles  for  direct  paths  are  under 
development.  A promising  scheme  obtains  an  initial  estimate  of 
the  depression  angle  from  a polynomial  fit  to  .sample  points  of 
the  Inverse  0(R) . These  estimates  are  refined  by  iterative 
techniques  until  the  error  is  within  satisfactory  bounds. 
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B.  PROPAGATION  LOSSES  FOR  THE  ECHO  LEVEL 

Once  the  above  scheme  Is  applied  to  find  both  the 
transmission  and  reception  ray  paths,  the  propagation  losses 
(including  attenuation,  spreading  losses,  and  losses  on  specular 
reflection  from  the  boundaries)  are  computed  as  set  forth  in 
Section  III. 

For  situations  in  which  ray  tracing  la  not  valid,  em- 
pirical formulae  will  be  implemented.  In  shallow  water,  for 
example,  ray  tracing  becomes  impractical;  the  COLOSSUS  formulae 
might  be  used  in  this  Instance.  Special  cases,  where  ray-tracing 
is  not  practical  are  handled  separately  and  are  not  discussed  in 
this  report. 

C.  DEVIATION  LOSS 

The  deviation  loss  t.s  a measure  of  tho  vimou\U  a signal 
is  diminished  by  the  army  pattern,  relative  to  signals  on  the 
beam  axis.  As  discussed  in  Section  IV,  the  pattern  Is  speci- 
fied by  two  angle  pairs  (H,  0),  and  - I'he  pair  wlil\  the 

subscript  "o”  refers  to  the  steering  (pointing)  direction,  and 
the  other,  to  the  observation  direction.  Vhe  ray  Inclination 
angle  Is  fixed  by  the  geometry  and  travel  lime  const  ritinls , 

For  recoption,  tho  ohservat  Ion  ar.lmvvlh  Is  spoctflod  bv 
the  angle  which  yields  the  mantmum  eclio  level,  denoleit  by 

(This  is  realistic,  since  It  yleld.s  t lie  most  probable  target 
bearing,  and  is  the  bearing  that  operational  systems  will  indi- 
cate.) For  cylindrical  sonars,  the  ohservat  Ion  a?,lmuth  Is  sim- 
ply  0^^  since  there  Is  no  beam  skewing.  For  C/F  sonars,  (unlike 
conventional  cylindrical  sonars)  the  most  probable  target  bear- 
ing depends  on  the  ray  depression  angle  , Itself  a function  of 
the  range.  (The  a/.lmuthal  angle  that  yields  the  maxlmvmi  echo 
level  is  given  by  0*  0^  (H^.,  0^.,  0,.^,),  with  H, 

and  0 held  constant.)  An  iterative  technique,  similar  to  the 

to 

one  described  in  Kefercnce  24  will  be  used  to  find  the  azimuthal 
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angle  which  maximizes  the  receiving  pattern  for  C/T  jonars . 

The  pattern  function  is  also  dependent  on  the  "point- 
ing” angles  (6^,  0^) . For  reception,  both  the  depression 
pointing  angle  0^^.!  and  the  azimuthal  pointing  angle  0^^,  are 
fixed  by  the  operator.  This  completes  the  specification  of 
the  reception  pattern.  The  transmission  pointing  angles 
(«of  also  are  decided  upon  in  advance.  For  transmission, 

the  observation  azimuth  i.s  determined  from  purely  geometric 
considerations . 

D.  THE  INSONIFICATION 

InsoniCication  is  defined  at  a point  "A”  as  the 
energy  flux  thruugh  a differential  area,  normal  to  the  ray 
path.  The  insonlf ication  process  can  be  characterized  by 
an  effective  pulse  length  and  a peak  intensity  for  each  direc- 
tion. This  pulse  length,  a constant  for  a fixed  depression 
angle,  is  chosen  so  that  the  prodvict  of  the  pulse  length  and  in 
tensity  function  describe  the  total  energy  flux  transmitted  in 
each  direction.  Propagation  losses  along  the  ray  path  from  the 
source  to  point  A are  taken  into  account  as  described  above. 

When  the  source  level  depends  upon  the  steering  direc 
tion,  the  source  inten.slty  will  be  represenied  by  a polynomial 
which  passes  through  sample  points  (c.g.  from  a cavitation- 
limited  source  level  computation,  or  experimental  data).  For  a 
single  transmitted  puise,  the  intensity  function  is  simply  the 
source  intensity,  modified  by  the  transmitting  puLtern.  When 
using  the  RDT  mode,  the  calculation  of  the  insonlf ication  is 
simplified  by  using  the  peak  intensity  at  the  point  "A”  and  an 
effective  pulse  length  to  yield  the  total  flux.  Since  the 
transmission  depression  steering  angle  is  fixed  by  choice,  the 
peak  source  intensity  for  a given  observation  direction  depends 
only  on  the  steering  azimuth.  An  iterative  teehnlque  locales 
the  azimuthal  beam  that  yields  the  peak  intensity,  allowing  the 
complete  specification  of  the  insonlf ication . 
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The  peak  intensity  function  (1^)  is  given  by; 


,(«t, 


«t)  - 


Max  Source 

Oot 


Intensity (9 


ot ' 


’ot) 


• (Transmitting  Pattern  ^ot*  ®t*  ^t^^] 


E.  THE  BISTATIC  TARGET  STRENGTH 

The  above  specifies  the  calculation  of  the  Insonlf- 
ication  and  propagation  losses.  Only  the  target  strength  for 
non-reciprocal  ray  paths  is  needed  to  complete  the  echo  level 
computation.  Information  on  the  target  strength  viewed  by  sonars 
using  a receiver  remote  from  the  source  Is  scanty.  Available 
sources  (See  references  21  and  22) , indicate  that  the  target 
return  is  quite  structured.  Averaging  over  small  solid  angles 
to  remove  this  fluctuation,  a mean  target  strength,  independent 
of  the  blstatic  angle  results  (geometry  is  illustrated  in  Figure 
5-1).  This  averaged  target  strength,  ranging  from  10  db  dead  on, 
to  20+  db  for  a broadside  view,  depends  only  on  the  aspect  angle. 
The  target  strength  represents  the  fraction  of  the  energy  inci- 
dent on  the  target  that  is  scattered  to  the  receiver. 

F.  THE  ECHO  LEVEL 

With  the  above  computations  completed,  the  echo  level 
is  given  by; 

ECHO  LEVEL  = SOURCE  LEVEL 

- DEVIATION  LOSS  (Transmission- 

- PROPAGATION  LOSS  (Transmission  path) 

- PROPAGATION  LOSS  (Reception  path) 

- DEVIATION  LOSS  (Reception) 

+ SIGNAL  PROCESSING  GAIN 

(All  quantities  are  expressed  in  db  relative  to  some  reference 
level  0) . 

The  echo  level  calculations  described  above  yield  the 
maximum  echo  level  that  the  system  will  obtain.  If  the  average 
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echo  level  Is  desired,  the  azimuthal  deviation  losses  must  be 
accounted  for;  these  depend  on  the  beam  spacings  used  in  trans* 
mission  and  reception.  (Note  that  the  beam  spacings  for  these 
two  modes  are  determined  Independently.  See  Reference  19  for 
further  discussion  concerning  the  selection  of  RDT  beam  spacings.) 


5-6 
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SECTION  VI 

BOUNDARY  REVERBERATION 

The  reverberation  from  the  horizontal  ocean  boundaries 
and  biological  scatterers  that  arrives  at  the  receiver  at  the 
time  of  the  echo  competes  with  the  echo.  Specular  reflections 
from  these  boundaries  have  discrete  travel  times  to  the  receiver, 
and,  generally,  are  not  likely  to  Interfere  with  echoes.  Non- 
specular  reflections  from  the  surface  and  bottom  were  character- 
ized In  Section  III.  The  discussion  of  the  deep-scattering 
layer  (DSL)  in  Section  III  Indicates  that  the  biological  scat- 
terers, at  the  frequencies  of  Interest  for  a C/P  sonar,  tend  to 
congregate  In  well-defined  layers.  For  the  purpose  of  estimating 
the  energy  scattered  by  the  DSL  to  the  receiver,  the  scatterers 
may  be  lumped  into  discrete  layers.  (This  corresponds  to  a 
numerical  integration  over  the  scattering  region  in  depth.) 

Accordingly,  the  three  forms  of  reverberation  may  be 
treated  by  a general  discussion  of  reverberation  from  a plane 
boundary.  The  contributing  area  on  the  surface  under  consider- 
ation Is  determined  by  simultaneous  application  of  temporal  and 
spatial  ccnstralnts.  With  the  contributing  area  defined,  the 
differential  expression  for  reverberation  is  integrated  to  yield 
the  reverberation  level  for  the  surface.  A sum  over  all  sur- 
faces (bottom,  surface,  and  DSL),  and  over  all  ray  paths  to  and 
from  these  surfaces  produces  the  total  reverberation  level.  In 
a discussion  of  boundary  reverberation,  it  is  convenient  to  put 
aside  a consideration  of  the  combinations  of  ray  paths  linking 
the  source,  reverberating  surface,  and  the  receiver,  and  concen- 
trate on  the  essence  of  boundary  reverberation.  Since  all  three 
types  of  boundary  reverberation  are  conceptually  the  same,  the 
following  discussion  of  boundary  reverberation  can  hold  for  all 
reverberating  surfaces,  for  any  combination  of  ray  paths  from 
the  source  to  the  reverberating  area,  and  from  the  reverberating 
area  to  the  receiver. 
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A.  SPECIFICATION  OF  THE  CONTRIBUTING  ANNULUS 

To  compete  with  the  echo,  the  energy  scattered  from 
the  reverberating  surface  must  arrive  at  the  time  of  the  echo 
(T_) . The  ray  path  associated  with  the  reverberation  must  start 
at  the  source,  strike  the  reverberating  surface,  and  end  at  the 
receiver.  Vector  notation  will  be  used  wherever  practical  in 
the  following,  to  point  out  that  the  results  are  Independent  of 
the  choice  of  the  coordinate  system.  Parameters  associated 
with  the  transmission  path  are  identified  by  the  subscript  "t"; 
those  related  to  the  reception  path,  by  a subscript  "r".  This 
leads  to  the  set  of  equations: 


where  s^  is  the  vector  between  the  source  and  receiver,  ('(^/r) 
is  the  vector  between  the  (source/receiver)  and  the  point  of 
contact  of  the  reverberating  ray  path  with  the  surface,  and  T 
denotes  travel  time. 

The  simultaneous  solution  of  these  equations  proi^eeds 
by  an  iterative  method.  (Both  travel  time  and  range  are  func- 
tions of  the  source  and  receiver  angles.)  The  Newton  iteration, 
operating  on  both  depression  angles,  is  convenient  once  an  ini- 
tial estimate  is  gotten  from  a solution  of  the  range  equation 
for  a point  on  the  reverberating  boundary  on  a line  with  the 
source  and  receiver.  An  Initial  estimate  of  the  depression 
angles  is  obtained  from  an  isovelocity  solution  chosen  to  give 
a total  travel  time  equal  to  the  echo  time.  In  general,  these 
ray  paths  will  not  meet;  holding  one  source  angle  fixed,  the 
other  is  varied  so  that  the  ray  paths  meet  on  the  reverberating 
surface.  This,  however,  yields  a travel  time  different  irom 
the  echo  time.  Then,  the  time  constraint  equation  is  modified 
to  cancel  the  difference  in  the  travel  times.  The  iteration  is 
continued  until  the  travel  time  and  range  error  (a  vector)  are 
within  satisfactory  limits. 
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An  Isoveloclty  ocean  restricts  the  reverberating 
area  to  an  annulus  having  the  shape  of  a cross-section  of 
an  ellipse.  The  ellipse,  with  the  source  and  receiver  as 
foci,  is  defined  by  the  echo  time.  The  wall  thickness  width 
of  the  annulus  is  determined  by  the  pulse  length.  (See  figure 
6-1)  In  the  usual  case  where  the  travel  times  of  interest 
greatly  exceed  the  pulse  length,  the  wall  of  reverberating  an- 
nulus is  quite  thin.  Since  the  insonification  is  nearly  con- 
stant across  the  narrow  width,  the  integration  of  the  differ- 
ential reverberation  may  be  reduced  to  a contour  integration. 
With  this  assumption,  one  of  the  integrals  of  the  surface  in- 
tegration is  done  implicitly.  The  surface  integral  reduces  to 
a line  integral  over  the  reverberating  contour.  Although  the 
results  which  follow  do  not  depend  on  the  reverberating  annulus 
having  the  shape  of  an  ellipse,  the  isoveloclty  geometry 
gives  a useful  Insight.  While  the  refractive  ocean 
will  change  the  shape  of  the  reverberating  annulus  somewhat, 
the  characteristics  of  the  reverberation  geometry  are  not  seri- 
ously altered. 

Once  a solution  to  the  range  and  travel  time  equations 
is  established  for  a point  on  the  line  connecting  source  and  re- 
ceiver (this  line  is  the  projection  of  onto  the  reverberating 
plane),  the  contributing  contour  is  traced  out  by  the  following 
scheme.  (The  technique  to  find  a point  on  the  connecting  line 
is  the  same  as  the  one  presently  used  in  TRG's  monostatic  cal- 
culations.) Figure  6-1  illustrates  the  geometry. 

1.  With  0^  and  0^  set  to  the  source  and  receiver  mutual 
relative  bearings,  the  resulting  2x2  set  of  equations 
for  horizontal  range  and  travel  time  is  solved  for  9^. 
and  0^,  as  for  the  monostatic  case,  by  the  iterative 
scheme  outlined  above. 
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FIGURE  6-1.  SIMPLIFIED  REPRESENTATION  OF  BISTATIC  REVERBERATION 
GEOMETRY  (THE  REVERBERATING  AREA  IS  SHADED) . 
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2.  The  transmission  depression  angle  (9^)  Is  Incremented  by 
a small  a9^.  The  time  constraint  gives  from: 


\9^  + 


nr  ■ 0 ' •'1 

r 


so  that 


3.  Given  the  depression  angles,  the  transmission  and  recep- 
tion azimuths  are  found  so  that  the  one-way  vector  ranges 
meet  at  a point  on  the  contributing  annulus.  Note  that  the 
reception  depression  angle  Is  chosen  for  9^  fixed  to  satisfy 
the  travel  time  constraint.  Since  the  azimuths  are  defined 
with  respect  to  the  sonar  cocrdiuate  system,  it  becomes 
necessary  to  resort  to  a coordinate  system  to  calculate  the 
azimuth  angles.  Ljt  denote  the  unit  vector  in  the  direc* 

tion  of  the  souxce/rece Lver , and  S be  the  unit  vector  in  the 
vertical.  Now,  |t  ,*  - §(^.2)  * vector  with  components 

only  in  the  horizontal  plane.  (Unsubscripted  vectors  stand 
for  both  transmission  and  reception.)  The  azimuthal  angles 
are  found  by  noting  that  cos  0 is  the  direction  cosine  of 
the  vector  with  the  heading  unit  vector. 

cos  0 . 

5.  At  any  point  on  the  reverberating  contour,  the  width  of 
the  annulus  is  found  from: 

AR  ,'(T)  - p(T“v)  , with  I standing  for  the  pulse 

length . 
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B.  THE  DIFFERENTIAL  REVERBERATION  INTENSITY 

The  reverbefating  area,  defined  Cor  a particular  set 
of  transtnlssion  and  reception  paths,  is  the  locus  of  points 
such  that  the  paths  from  the  source  to  the  surface  and  from  the 
surface  to  the  receiver  have  a travel  time  equal  to  the  echo 
time (within  a few  pulse  lengths).  The  differential  reverber- 
ation calculation  is  conceptually  the  same  a.s  the  echo  calcul- 
ation, when  the  target  is  taken  to  bo  a differential  area  on  the 
reverberating  boundary.  The  .scattering  coefficient  |i,  a func- 
tion of  the  direction  cosines  with  re.spect  to  the  normal  to  the 
surface  is  chosen  to  conform  to  observed  data  lor  the  surface 
of  interest. 

Then,  the  general  expression  for  the  differential  re- 
verberation intensity  due  to  a reverberating  boundary  is; 

dl^,  n(e[.  0)V,.(9,.,  U)dA 

where 


M is  the  scattering  coelflcleni  per  vtnll  area, 
is  a grayling  angle  at  the  Innnuiary 

V is  a / limit  ha  1 angle  between  ihe  Incident  and 
scattering  direct  tons 

K Is  the  ray  bundle  spreaillng  loss  lactor 
Vj.  Is  the  receiving  Inti'nsltv  tvinctlon 

I,  L.s  the  transmitting  source  Inlensltv  lunction, 
which  lor  a single  pulse  Is  given  by 


wU.’i'i'  (9,,,,  0^,,)  Is  ru.’.l, 

and  lor  RDT  is  given  lyv 

where  lixed. 
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where  : 

is  the  transmitting  intensity  pattern  function, 

is  the  source  intensity  function,  and 
o 

dA  is  the  differential  area  of  concern. 

Tlie  selection  of  a coordinate  system,  useful  for  the 
numerical  Integration  of  the  rever'  i-.-ation,  depends  \ipon  the 
character  of  the  insonif ication  of  the  reverberating  area.  This 
is  under  investigation. 

C.  THE  NUMERICAL  INTEGRATION  FOR  THE  REVERBERATION 

Since  a closed-form  expression  for  the  differential 
reverberation  intensity  has  not  been  found  for  the  general  case, 
a numerical  integration  over  the  contour  is  used  to  calculate 
the  reverberation.  Previous  experience  with  numerical  integra- 
tions of  this  nature  indicates  that  adaptive  schemes  are  not 
very  helpful;  they  seem  to  require  a number  of  sample  points 
approximately  equal  to  twice  the  number  of  peaks  in  the  inson- 
Ification  of  the  annulus.  A straight-forward  scheme,  based  on 
the  upper  and  lower  sums  that  bound  the  Riemann  integral,  appears 
promising,  Tills  scheme  has  the  advantage  that  a truncation 
error  estimate  is  available  from  a comparison  of  the  upper  and 
lower  sums . 

Upon  examination  of  the  bistatic  insonificat  Ion  ol  tlic 
reverberating  areas,  an  analytic  approximation,  patterned  after 
the  one  detailed  for  the  monostatic  case  in  Appendix  A of  Refer- 
ence 24  may  be  developed  to  replace  the  numerical  integration 
wherever  possible. 

D.  ECHO -TO -BACKGROUND  RATIO 

With  the  calcultalons  presented  above  completed,  the 
echo-to-background  ratio  is  easily  computed  from: 

E/B  - Echo  Level  - (Total  Reverberation  ! Sc  If -Noise) . The 
estimation  of  self-noise  is  de.scribed  in  Section  IV,  the  echo- 
level,  in  Section  V,  and  the  reverberation  (summed  over  all  types 
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and  orders),  in  this  section.  The  noise  and  reverberation  compo- 
nents are  assumed  to  have  random  phases  and  their  powers  added 
to  obtain  total  background  level.  These  echo-to-background 
ratio  calculations  provide  the  means  for  a study  of  blstatlc 
echo  ranging.  A full  understanding  of  the  nature  of  echo 
ranging  with  systems  using  a remote  receiver  may  lead  to  more 
fruitful  employment  of  present  sonar  systems,  and  provide  a 
guide  for  the  design  of  future  ones. 
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